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ABSTRACT 



o. 

, We aim to construct a general relativistic radiative transfer formulation, applicable to particles with or without mass in astrophysical 

settings, wherein ray-tracing calculations can be performed for arbitrary geodesies for a given space-time geometry The relativistic 
radiative transfer formulation is derived from first principles: conserving particle number and phase- space density The formulation 
is covariant, and transfer calculations are conducted along particle geodesies connecting the emitters and the observer The geodesies 
are determined through the space-time metric, which is specified beforehand. Absorption and emission in the radiative transfer cal- 
culations are treated explicitly The particle-medium interaction is evaluated in the local inertial frame, co-moving with the medium. 
Relativistic, geometrical and optical depth effects are treated self-consistently within an integral covariant framework. We present a 
self-consistent general relativistic radiative transfer formulation with explicit treatment of emission and absorption. The formulation 
is general and is applicable to both particles with mass and without mass. The presence of particles has two major effects: firstly 
the particle bundle ray is no longer along the null geodesic, and secondly the intensity variation along the particle bundle ray is 
reduced by an aberration factor. The radiative transfer formulation can handle 3D geometrical settings and structured objects with 
variations and gradients in the optical depths across the objects and along the line-of- sight. Such scenarios are applicable in calcula- 
tions of photon emission from complex structured accretion flows around black holes and neutrino emission from remnant neutron 
tori in neutron-star mergers. We apply the formulation and demonstrate radiation transfer calculations for emission from accretion 
tori around rotating black holes. We consider two cases: idealised optically thick tori that have a sharply defined emission boundary 
surface, and structured tori that allow variations in the absorption coefficient and emissivity within the tori. We show intensity images 
0^ | and emission spectra of the tori obtained in our calculations. Our findings in the radiative transfer calculations are summarised as 

follows, (i) Geometrical effects, such as lensing-induced self-occulation and multiple-image contribution are much more significant 
in accretion tori than geometrically thin accretion disks, (ii) Optically thin accretion tori show emission line profiles distinguishable 
from the profiles of lines from optically thick accretion tori and lines from optically thick geometrically thin accretion disks, (iii) The 
line profiles of the optically thin accretion tori have a weaker dependence on the viewing inclination angle than those of the optically 
thick accretion tori or accretion disks, especially at high viewing inclination angles, (iv) Limb effects are present in accretion tori 
with finite optical depths, due to density and temperature stratification within the tori. We note that in accretion flows onto relativistic 
compact objects, gravitationally induced line resonance can occur. This resonance occurs easily in 3D flows, but not in 2D flows, such 
as a thin accretion disk around a black hole. 

Key words, accretion, accretion disks - black hole physics - galaxies: active - line: profiles - radiative transfer - relativity 
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1 . Introduction to emission from the two parts of the disk, which have opposite 

! • ' , ^ . . t , . . , t . / . projected line-of- sight velocities. Double-peaked optical lines 

. £ i T1 \ e * " r f fission observed in active galactic nuclei (AGN) haye been obseryed i n a varie tv f binary systems (e.g . black 
^ and black hole binaries is be lieved to be powered by accretion ho i e X-ray binar ies. [Johnston et all flM iMarsh et all Il994t 
K£ . of material onto bl^k hol es igalBSldlpq iLynden-Belfl gp | S oria et alJ[T999t |Wu et al. 200lDouble- peaked optical lines 
; I Shakura&Sunyae^ BlIl. The accreting hot plasmas rotating ^ al so seen in a small fra ction of^NjPuchnarewic^iriD 



around the black hole form a disk or a torus. It has been sug- rj^. | Er acleous & Haloed l200l IStrateva et alj £o06h . These 

gested that dense remnant neutron ton can also be formed around double . peaked lines can be explained in a Newtonian fram ework 

compact objects which could be a very massive n eutron star or as described in Sma g ^ (see also | H orne & Marshi ri986i. 

ajjke k hole after tvvo neutron stars merge (see | Shibata et al. | Double . peaked lines have also been observed in the X - ray S pec- 

Ml Baiotti et al. M| Rezzolla et alM 3). ^ such systems, ^ of ^ ^ Broad a tric double ^ e aked 

the influence of curved space-time is significant. It affects the ra- Fe Ka ^ we re found in the spec tra of a number of AGN (e.g. 

diative transport of particles in the ac cretion flow as wel as th e MCG . 6 . 3 o. 15 . ffa n ^7 t al . 1995) . The X -rays of AGN are be- 

hydrodynamics of the flow itself (see |Novrko y & Thorne| [l9g). lieyed tQ Qri inate from ions dose tQ me central black 

Emission from accretion disks around compact objects has ^ whefe ^ accretion flow Js M M relativistic and the 

been investigated for several decades now. Emission lines from { h The emissions from different rts of the accre . 

geometrically thin accretion disks aroun d gravitating objects are tion flow ^ therefore boQSted differentiall Various relativistic 

expected to have two peak s <lSma!dll969|). The peaks correspond eflfects alsQ cause additional differ ential broadening and distor- 

* Send offprint requests to: Z. Younsi ti° n °f an y ^ me emission. As such, the emission lines from the 
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inner regions of relativistic accretion disks around black holes 
have a very broad pr ofile, with an extended red wing and a pro- 
nounced blue pea k (ICunninghaml 119751: iReynolds et al.l 119991: 
iFabian et"ai1l2000l) . Black holes with faster spins would give rise 
to relativistic accretion lines with a broader red wing, because 
the inner boundary of an accretion disk around a maximally ro- 
tating black hole can extend very close to the black hole event 
horizon. 

The emission from accreting gas and outflows in the 
vicinity of black holes is subject to Doppler shifts, lensing, 
gravitational time-dilation, and other relativistic dynamical 
effects. There have been numerous calculations of relativistic 
(phot on) lines from accre ti on disks and tor i around black holes 
(e.g. I Cmmin^haml [1 9751: iGerbal & Pejajl 119811: [Fabian et al.l 

iFanton et al. l 
iFuerst & Wul 




Stella! ll99Ct iKoiimal [l99lt iBao 



Reynolds et al 



1999; Fabian et al 



Beckwith & Donel l2QQ4al 120051: ICadez & Calvanil 



Schnittman & Rezzolla~|SQl~ Fuerst & Wu 



Dexter & Agol 2009; Soch ora et aD 1201 ll Ivhicent et al.l 



20071: 



Wang & Lil 12012b . The three most common methods for cal- 
culating relat ivistic line profiles are (i) the transfer function 



method (e.g. ICunninghaml 119751: IFabian et al l 120 00). (ii) the 
elliptic function method (e.g. iDexter & AgoTl2009l) and (iii) 
the dir ect geodesic integratio n metho d (e.g. IFuerst & Wul 
20041: ISchnittman & Rezzollal 120061: IFuerst & Wul 120071: 
Anderson et al .1120 lot I Vincent et al.ll201ll) . The transfer function 



method and the elliptic function method are efficient for the 
calculation of emission from thin axisymmetric optically thick 
accretion disks but are not applicable to systems that lack 
the appropriate geometry and symmetry. The direct geodesic 
integration method is a brute force approach, and less restrictive 
in this context, compared to the other two methods. It works 
well with any three-dimensional (3D) accretion flow, e.g. 
time-dependent accretion flows from numerical relativistic hy- 
drodynamic simulations, and can also handle opacity variations 
within the system. In most of these relativistic calculations, the 
focus was on the investigation of line broadening due to rela- 
tivistic effects. While relativistic ray-tracing of photons in strong 
gravity and the corresponding calculations of emission line pro- 
file broadening have been investigated in various astrophysical 
settings for decades, there have been only a few studies, often in 
restricted settings, of the opacity effects due to self- absorption, 
emission and scatteri ng within the accretion flows and along 
the line of s ight (e.g. IZane et al.l[l996L iFuerstl 120061: IWu et al.l 
120061 2008; D olence et al.l 120091) . Covariant radiative transfer 
calculations of emission from accretion flows in more general 
settings with an explicit treatment of absorption, emission and 
scattering are lacking. There is also no corresponding covariant 
radiative transport formulation applicable to both relativistic 
particles without mass (photons) and with mass (e.g. neutrinos), 
in a general astrophysical setting in the present literature. (Note 
that particles with mass do not follow null geodesies, and the 
relativistic formulation and the corresponding ray-tracing need 
to be modified.) 

In this work we present the general covariant formulation for 
radiative transport of relativistic particles. The radiative trans- 
fer equation is derived from the Lorentz-invariant form of the 
conservation law. It explicitly includes emission and absorption 
processes. The formulation is not restricted to general relativis- 
tic radiative transfer of photons. It is general and can be ap- 
plied to both relativistic particles without mass (e.g. photons) 
and with mass (e.g. relativistic electrons and neutrinos). The ra- 
diative transfer recovers its conventional form in the Newtonian 
limit. We carried out demonstrative radiative transfer calcula- 



tions of emission from model accretion tori around rotating 
black holes with different optical thicknesses. Our calculations 
show the convolution of geometrical and optical depth effects, 
such as limb darkening and brightening caused by optical-depth 
and emissivity variations, multiple image contribution, and ab- 
sorption and self-occultation induced by gravitational lensing, 
which are characteristic of 3D structured flows in strong gravity 
environments. The paper is organised as follows. In Sec. 2 the 
covariant radiative transfer equation is derived and its solution 
discussed. In Sec. 3 we show the construction of two torus mod- 
els, one with a specific sharp emission boundary surface, and an- 
other with a stratified density and temperature structure. In Sec. 
4 we present the results of radiative transfer calculations for the 
torus models. We briefly discuss the astrophysical implications 
of our calculations and the possible occurrence of gravitationally 
induced line resonance in 3D accretion flows near black holes. 



2. Covariant radiative transfer formulation 

The emission and absorption processes may be considered as 
sources and sinks in the medium through which the ray bundle of 
particles is transported. In their absence, the number of particles 
along a ray bundle is conserved. In the presence of emission and 
absorption, the radiative transfer equation can be expressed as 



r ~ CXyly + Jy , 

as 



(D 



where I v = I v (s) is the specific intensity of the ray at a frequency 
y, and a v and j v are, respectively, the absorption and emission 
coefficients at a frequency v. By introducing the variable 



I a v (s') ds' 

J s 



(2) 



which is the optical depth (optical thickness of the medium be- 
tween s and so), we may rewrite the radiative transfer equation 
as 



^ = -I v + ± = -I v + S v 

ciTy a v 



(3) 



where S v = j v /a v is the source function. Direct integration of 
the equation yields 

I v (s) = I v (s )e- T *+ P j v (s')e- ( ™- Ms ' }) ds' 

J S 

= I v (0) e" T * + f V S v (t' v ) e- (T ^W v , (4) 
Jo 

where the constant I v (sq) (= / v (0)) is the initial value of the spe- 
cific intensity. While optical depth, which is a scalar quantity, 
is invariant under Lorentz transformations, the radiative transfer 
equations in the conventional form (equations Q and 0) are 
not. 

We now show that a covariant formulation for radiative trans- 
fer can be derived from the conservation of phase space vol- 
ume and the conservation of particle number. Construct a phase 
space volume <V threaded with a small bundle of particles. In 
the co-moving frame, these particles occupy a spatial volume 
element d 3 x = dx dy dz and a momentum volume element 
d 3 p = dp x dp y dp z . Liouville's theorem ensures that the phase 
space volume, given by d*V = d 3 xd 3 p, is unchanged along the 
affine parameter A, i.e. 



~dX 







(5) 
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(Mis ner et al .11 19731) . This together with the conservation of the 
number of particles dN within the phase space volume element 
implies that the phase space density 



f{x\p { ) 



dN 



(6) 



is invariant along A. 

For relativistic particles, \p\ = E, and d 3 p = E 2 dEd£l (here 
and hereafter, unless otherwise stated, we adopt the convention 
that the speed of light in vacuum c - 1). The volume element 
of a bundle of relativistic particles is also given by d 3 x = dA dt, 
where dA is the area element of the bundle. Thus, we may ex- 
press the phase space density of a bundle of relativistic particles 
as 



dN 



f( X ' p) E 2 dAdtdEd£l ' 
The specific intensity of a ray (bundle of photons) is simply 
EdN 



(7) 



(8) 



(9) 



dAdtdEdQ 

fsee lRvbicki & Lightmanll 19861) . It follows that 

T= k = II 

y 3 £ 3 

is a Lorentz invariant quantity. We denote I as the Lorentz- 
invariant intensity, and it can also be regarded as the occupa- 
tion number of particles in the phase space for a particle bun- 
dle. We may also obtain the corresponding Lorentz-invariant ab- 
sorption coefficient^ = va v and the Lorentz-invariant emission 
coefficient r] = j v /v 2 . These two coefficients, as seen by the ob- 
server, are related to their counterparts in the local rest frame of 
the medium via va v = vo#o,v and j v /v 2 = Jo,v/Vq respectively, 
where the subscript "0" denotes variables measured in the local 
rest frame. 



2.1. Photon and relativistic massless particle 

For photons (or a massless relativistic particle) k a k a = 0, where 
k a is the (covariant) 4-momentum. Consider a photon propagat- 
ing in a fluid with 4- velocity . The photon's velocity in the 
co-moving frame of fluid, v^, can be obtained by projecting the 
photon's 4-momentum into the fluid frame, i.e. 

/ = P a ?k a 

= kP + (k a u a )if , (10) 

where F*P = g a P + u a u^ is the projection tensor and g a P the 
space-time metric tensor. The variation in the path length s with 
respect to the affine parameter A is then 

% = -'^"L 

a A Mobs 



= - ^gap(k a + (k p uP)u a ){kP + (k a u a )uP 
= - JkpkP + (k a ii a ) 2 upuP + 2(k a u a ) 2 \ 



^obs 



— k a u 



Mobs 



(11) 



Here we have used [-,+,+,+] convention for the signature of 
the space-time metric. For a stationary observer located at infin- 
ity, ppu 13 = -E b s . The relative energy shift of the photon be- 
tween the observer's frame and the comoving frame is therefore 



-i _ vo _ -k a u a \ x _ k a u a \ x 
v E obs huP\x obs 



(12) 



Making use of the Lorentz-invariant properties of the vari- 
ables J, x an d T], and of the optical depth r v , we may rewrite the 
radiative transfer equation in the following form 



^ = -I + «- = -I + S, 

dTy X 



(13) 



where S = rj/x is the Lorentz-invariant source function. All 
quantities in equation (TT3l) are Lorentz invariant, and hence the 
equation is covariant. As dr v = a v ds, we have 



— = -a v I + — . 
ds v 5 

It follows that 

dI r ai I t , hv \ 

dX = ~ kaU u \~ a ^ I+ ~ I ' 



(14) 



(15) 



where oro, v = ^o(^> y ) an d jo, v = joi^ 3 , v), where, as before, "0" 
denotes variables that are evaluated in the local rest frame, for a 
given y. In ray-tracing calculations we specify an observer fre- 
quency v, determine all required variables at this location, and 
find how these variables change in different reference frames 
through the radiative transfer equation. The solution to equation 



1(A) = J(^ )e 



TAX) 



(cf. iBaschek et al.lll997l: iFuerst & Wull2004l) . where the optical 
depth is 



T V (A) = - I a^ v (A)k a u a \ x >dA' . 

In terms of the optical depth, 

J(r v ) = /(To)e-^ + r 5«)e- (Ty -^ dr' v 

Jt 



(17) 



(18) 



For a distant observer, -k a u a \x obs — > E, the observed photon 
energy, which may be normalised to unity. The radiative transfer 
equation can be expressed as two decoupled differential equa- 
tions 



dr v _i 
~dA =r " 



dl_ 
~dA 



i J0,v\- T 



(19) 
(20) 



These two equations are more useful in practical relativistic 
radiative transfer calculations because they allow the efficient 
computation, through a simple Eulerian method, of the optical 
depth along a ray, regardless of whether the ray-tracing is exe- 
cuted forward or backwards in time. This is then used to compute 
the intensity along the ray. 

2.2. Particles with mass 

For massless particles, contraction of the 4-momentum gives 
k a k? = 0, but for particles with a non-zero mass m, it gives 
PaP a = -m 2 . The presence of the particle mass modifies the 
particle's eq uations of motion, changing the geodesi es from null 
to time-like (ICartedll968l: iBover & Lindquistll 19671) . Moreover, 
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a covariant particle flux is mass dependent. The radiative trans- 
fer equations derived for massless particles (equations (TPfl) and 
(fT5l) ) are therefore not applicable for particles with mass, such as 
neutrinos or relativistic electrons. Nevertheless, the formulation 
obtained for massless particles can easily be modified to include 
the effects due to the particle mass. We follow a similar proce- 
dure as in the case of massless particles in our derivation. We 
express the variation in the path length with respect to the affine 
parameter as 

ds / 1 

77 = -J8a/3(p a + (p^)u a )(pP + (p a u a )ufi)\ , (21) 

CLA " \Aobs 

analogous to equation (ITTI) in the case of massless particles. 
Because p a p a + 0, we have 



els 
~dX 



■^P(3pP + {PaU a ) 2 



Aobs 



~aJ(Pi 



(22) 



Insert this into equation ([T3l) , which does not depend on the par- 
ticle mass explicitly. After some algebra, we obtain the general 
covariant transfer equation for relativistic particles 



dl 
dA 



\ 



( \ 


2 




m 


p a u a 


J 


VruP\ . 

V vp \Aobs ' 







70,1 



(23) 



For a stationary observer located at infinity, ppuP = -E. 
Equation (l23t differs from equation (1131) by an aberration factor 
yjl - (m/E) 2 . This factor reduces the intensity gradient along 
the ray. Because it approaches unity when m — > 0, the radiative 
transfer equation (1131) is the radiative transfer equation (1231) in 
the zero-mass limit. 

3. Accretion tori around black holes 



We can see in equation (1231) that the essential components of the 
covariant radiative transfer formulation are the emission coeffi- 
cient, the absorption coefficient, the relative energy shift of the 
particles with respect to the medium, and the aberration factor, 
which are all evaluated along the particle geodesies. The particle 
geodesic is obtained by solving the equation of motion of parti- 
cles for the specified space-time metric with the location of the 
observer assigned. The local values of the emission coefficient 
and the absorption coefficient can be calculated when the ther- 
modynamic conditions of the medium are given, and the relative 
energy shift and the aberration factor can be determined when 
the hydrodynamic properties of the medium are known. 

Here we demonstrate the application of the covariant radia- 
tive transfer formulation and calculate the emission from 3D ob- 
jects in a gravitational field. We consider accretion tori around 
rotating black holes. The tori may have several emission compo- 
nents with different optical depths. Accretion tori are 3D objects 
with internal structure. They are in contrast to optically thick, 
geometrically thin accretion disks, which are 2D objects where 
explicit covariant radiative transfer is unnecessary in determin- 
ing how the radiation propagates and is modified inside the disk. 

3.1. Modelling accretion tori 

For accreting objects, the accretion luminosity L acc roughly 
scales with the mass accretion rate M as L acc — sM. The con- 
version parameter s ~ 10 20 erg g" 1 for black holes. The forma- 
tion of geometrically thin accretion disks around a black hole 



requires that the radiation pressure force in the accretion flow 
is much smaller than the local gravitational force exerted by the 
black hole. This condition is usually satisfied when M is suf- 
ficiently low, such that L acc is much lower than the Eddington 
luminosity, which is given by L E dd = 1.4 x 1O 38 (M/M ) erg s" 1 , 
where M is the mass of the accretor, and M the solar mass. 
Because L acc increases with M, high M implies high L acc and 
hence a large radiative pressure within the disk where the radia- 
tion is liberated. When L acc approaches L E dd, the radiation pres- 
sure force become comparable to the local gravitational force, 
and the accretion disk thereby inflates and become a torus (see 
e.g. lFrank etaDl2QQ2h . 

In general, full knowledge of the fluid viscosity is required 
in determining the structure and hydrodynamics of the accre- 
tion torus. However, in accretion tori the angular momentum 
transport is non-local. The process cannot be parametrised with 
a local viscosity, as in the case of modelling geometrically 
thin accre tion disks, where the a-v i scosity prescription is of - 
ten used dShakura & Sunvaevl Il973l: lAbramowicz et alJll988h , 
It is believed that the angular momentum transport in accretion 
disks/tori is mediated by tangled magnetic fields that permeate 
the flow (e.g. magn eto-rotational instability (MRI), see lHawleyl 
l2000l: iBalbus 2003). In principle this magnetic viscosity and the 
flow hydrodynamics need to be determined simultaneously and 
self-consistently. Nevertheless, certain phenomenological pre- 
scriptions are proposed to bypass the viscosity calculations, e.g. 
assuming an angular momentum distribution within the torus in- 
stead of solving for the distribution. With this, the structure of 
the accretion torus can be determined by solving only the re- 
maining hydrodynamic equations and the equation of state. For 
the purposes of this work we considered this phenomenologi- 
cal approach and constructed accretion tori assum ing a specific 
angular velo city profile within the to rus (see also IFuerst & Wul 
12004 120071) and kbramowiczl (12005b . With the angular veloc- 
ity profile specified, we determined the density and entire flow 
profiles in terms of certain normalised variables. The resulting 
accretion torus was then rescaled, using the results fro m accre- 
tion to r i/disks obtaine d by numerical MRI simulations (lHawleyl 
(l2QQQh : iBalbusI (120031) ). The torus model constructed as such is 
able to capture the geometrical aspects of the MRI accretion tori 
and the physical conditions within the accretion flow. 

3.2. Emission surface of rotationally supported torus 

To calculate the emission from an opaque accretion torus, we 
need only to specify the torus' boundary surface and its phys- 
ical conditions. The simplest model that allows us to specify a 
boundary emission surface is a rotationally supported torus. In 
it the total pressure force is balanced by the 4-acceleration in an 
arbitrary fluid element. The inner boundary of the torus is de- 
fined where this balance breaks down. With an appropriate pa- 
rameterisation of the angular velocity profile (Q, as a function of 
position in the torus) we can derive the 4-acceleration and hence 
obtain the pressure force. Tracing the isobars gives the isobaric 
surfaces in the torus. The process is essential in constructing the 
gradient contours of the 4-acceleration in the r - 6 plane (in the 
(t, r, 6, cp) spherical co-ordinate system). The torus boundary sur- 
face is simply the outermost allowed isobaric surface. 

We consider a stationary, axisymmetric, rotationally sup- 
ported accretion torus. The symmetry axis of the torus is aligned 
with the spin vector of the central black hole. The 4-acceleration 
of the flow is derived from the 4-velocity 



(24) 
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Fig. 1. Cross-sections of the boundary surfaces of model rota- 
tionally supported tori in cylindrical coordinates with z = r cos 6 
(the equatorial plane of the tori is z = 0). The top panel shows the 
tori around Kerr black holes with spin parameters a = 0, 0.5 and 
0.988. The angular velocity profile index of the tori n = 0.21, 
and the Kelperian radius - 12 r g . The bottom panel shows the 
tori with angular velocity indices n = 0.14, 0.17 and 0.232. The 
Keplerian radius of the tori is /> = 8, and the black-hole spin 
parameters = 0.998. 



where Pj^ is the Christoffel symbol. Because the poloidal com- 
ponents of the flow are dynamically unimportant, we set the 4- 
velocity as u a = (V, 0, 0, u*). It follows that u a ^ vP = and 



(25) 



In Boyer-Lindquist co-ordinates, the r and components of the 
4-acceleration are 



A 



/X-2r 2 \/ \ 2 
Ml — Ui-asm 2 6<pj + rsin 2 (9 2 



(26) 



a e = "La 9 = - sin 2(9 \^[ai - (r 2 + a 2 )<pf + ^0 2 j , (27) 

where A = r 2 - 2Mr + a 2 , £ = r 2 + a 2 cos 2 6, and M and a are 
respectively the mass and spin parameters of the central gravitat- 
ing body. Setting a a u a = yields a set of differential equations 
for the isobaric surfaces. It is more convenient to express the 
equations as 



dr 



if/ 2 



dO _ -if/] 

* ^ 2 +A < A 2 ' 



(28) 
(29) 



for the calculation and construction of the torus boundary sur- 
face. Here £ is an auxiliary variable, and the functions if/\ and if/2 



are given by 
<Ai = Ml 



I? 



j^Q" 1 -asinoj 



+ r sin 



(30) 
(31) 



<A 2 = sin2^|^[aQ- 1 -(r 2 +a 2 )] 2 + ^J. 

Solving equations ([28b and ([29b . with the inner boundary radius 
of the torus and the specified angular velocity profile gives the 
torus boundary surface. 

We consider an angular velocity profile with the form 



Q(r sin 0) = 



(32) 



(r sin#) 3 / 2 + aVMWin6V 

following iFuerst & Will (l2004l l2007h . where r k is the radius 
on the equatorial plane at which the material circulates with a 
Keplerian velocity. The differential rotational velocity gradient 
gives rise to an implicit pressure force, supporting the torus ma- 
terial above and below the equatorial plane. The parametrisation 
using the variable r sin ensures the constant density and pres- 
sure surfaces coincide in the Newtonian limit and a polytropic 
equation of state is applicable for the flow. The index parame- 
ter n is crucial for regulating the pressure forces, thus adjusting 
the torus' geometrical aspect ratio. Its property is similar to that 
of the q index of the von Z eipel parameter in the study of sta- 
bility of accretion disks (seelChakrabarti 19851 iBlaes & Hawleyl 
Il988l) . Generally, n « q - 1.5, with the relation being exact for 
Schwarzschild black holes. Tori with q > V3 are unstable. In 
the Newtonian limit tori with n = 0.232 are marginally stable. 

The angular velocity and angular momentum of the flow are 
Q = uf/u* and I = -u^/u t respectively, and the redshift factor is 
given by 

A = u* = [-{g tt + mgt* + &g H )]~ l/2 , (33) 
and the energy per unit inertial mass of the flow material is 



U 



-u t 




(34) 



(35) 



+ l 2 gu ' 

These two quantities are related via 

Zero values for the denominators in equations (f33l) and (f34l) cor- 
respond to the locations (the photon surface) where the local 
flow speeds reach the speed of light. 

In this model we need to specify the condition for the inner 
boundary radius. We take it as the intersection of the isobaric 
surface with either the orbits of marginal stability or the limiting 
surface of photon orbits, whichever has a larger radius. Usually 
the photon surface is within the marginally stable orbit. The in- 
ner boundary of the torus is therefore in general determined by 
the outermost surface of the torus, which satisfies dU/dr = 0. 
Solving this gives 



2aMsin 4 6> 



r 2 +a 2 + 



a 2 Mr sm z 6\L-2r 2 



+ sin 2 6 



6aM 2 rsin 2 6 



+A sin 2 6>Q— 
or 



M\ 6Mr(r2+a2) + 3A- 



S 

E-2r 2 
-Mil - 



X 2 
■2r 2 



Z 2 



+ r 1 



2Mr\Z 

~ r 



■2r 2 



Z 2 



= 0. 



(36) 
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which when solved for r, yields the inner edge of the torus (for 
a particular choice of t> and n). Equations (f28l) and (|29t for the 
torus surface are now readily integrated. 

Figure 1 shows the boundary surfaces of rotationally sup- 
ported tori with various system parameters. The shape of the tori 
is determined by the rotational velocity index of the torus n and 
the black-hole spin parameter a. When is fixed, the vertical 
thickness of the torus increases with a but decreases with n. The 
degree of self-occulation of an optically thick torus and hence 
the spectral properties of the emission depend on the aspect ra- 
tio of the torus and the viewing inclination. Note that the tori in 
Figure 1 are purely rotation- supported. The thermal pressure of 
the torus gas and the radiative pressure of the emission from the 
gas have not been included in their construction. The presence 
of gas pressure and radiation pressure will modify the aspect ra- 
tio of the tori. In the next subsection we will consider a more 
general situation, which includes the gas pressure and radiative 
pressure, and construct tori with internal density and temperature 
structures. 

3.3. Pressure supported torus structure 

Accretion tori resemble stars that have an atmosphere with an 
optical depth gradient. While emission from an opaque torus 
stems from an unobscured skin (surface) layer of the torus, emis- 
sion from a translucent or an optically thin accretion torus is 
formed by the emission contribution from all regions within the 
torus. The thermodynamic and hydrodynamic structures of the 
torus determine the spectral properties of its emission, and so 
they must be determined prior to performing the radiative trans- 
fer calculations. 

To model the internal struc ture of the accreti o n tori , we 
adopted a prescription given by lAbramowicz et al. I (fl978h and 
iKozlowski et al.l (Il978l) . We considered the tori as stationary and 
axisymmetric. They consist of a perfect fluid, and the stress- 
energy-momentum tensor of the flow is given by 



T ap = (p + P + e)u a u p + Pf , 



(37) 



where P is the pressure, p is the density, and e is the fluid internal 
energy. Because T a ^.^ = 0, we have 

(p + P + e)fiU a vP + (p + P + e){u a fi vP + u a if. p ) + P#g afi = .(38) 

Projecting perpendicular to the velocity with the projection ten- 
sor P a P yields the momentum equation 



(p + P + e)u a y + P#fP = . 



(39) 



Because the torus is stationary and axisymmetric, it has neg- 
ligible poloidal flow components. Hence, 



Differentiating u a u a = -1 gives 

(u a U a ), S = g a ^5U a U^ + 2u a U a , S = . 

Thus, 



u u a s — ~ ~u 



SaB,S 



It follows that 

Ma-fill^ = uPup, a = lid a U t + U^daUfj) , 



(40) 



(41) 



(42) 



(43) 



where d a is the gradient in the x a direction. Recall that Q = uP/u* 
and / = -Uf/ut. Therefore, we have 



The gradient of / is 



Q 



1 -zn 



uf U t 
It follows that 

— = —o a u t + u d a u t + wd a u<p . 

1 - III u t 

Equation (l43t can now be expressed as 

a OdJ 1 

Uq-rW — — o a u t . 



Thus, 



d a P 



p + P + € 



i-ia u t 



= d a \n{u t ) - 



1 - ZQ 



(44) 



(45) 



(46) 



(47) 



(48) 



(cf. equation (7) in lAbramowicz et al. I (fl978h for the accretion 
torus). 

An equation of state is required to close the system of equa- 
tions and the gas within the torus is assumed to be barotropic. 
Owing to the complexity of the Kerr metric, the solution must 
be computed numerically. However, in the special case of Z = 
constant, corresponding to a marginally stable torus, there is an 
analytic solution 



f 

Jo 



dP' 



p + P' + 6 



= ln(M,) - ]n(ut)i, 



(49) 



where \n(u t )i nner is evaluated at the inner edge of the torus. 

The transparency of the torus to radiation requires that the 
local emissivity and opacity within the torus must be specified 
explicitly. As such, the velocity, temperature and density struc- 
ture of the torus must be determined prior to the radiative transfer 
calculations. 

The total pressure within the torus is the sum of the gas pres- 
sure and the radiation pressure, i.e. P = P gas + P m d-> where 



gas 



pkT 

pm H 
Act 



■PP., 



Prad = —T 4 = (l-fS)P. 

3c 



(50) 
(51) 



Here k is the Boltzmann constant, p the mean molecular weight, 
run the mass of hydrogen,/? the ratio of gas pressure to total pres- 
sure, and cr = n 2 k 4 /60fi 3 c 2 is the black-body emittance constant. 
Eliminating kT in the above equations yields 



P = he 



45(1-/?) 



7i\pm H pY 



1/3 



(52) 



For a poly tropic equation of state P = /cp r , the internal energy 
is related to the pressure by e = P/T - 1 . Equation (l52l) implies 
T = 4/3, and k = ftc[45(l - J3)/7r 2 (pm H /3) 4 ] 1/3 . 

Combining the momentum equation (|39t for a perfect fluid 
with the polytropic equation of state gives 



P+^—jP)a a = -P^ 



„O0 



(53) 
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Fig. 2. Cross-sections show the density, pressure and tempera- 
ture contours in model translucent tori (panels from top to bot- 
tom). The angular velocity profile index of the tori n = 0.21, 
the Keplarian radius = 12 r g . The black-hole spin paramter 
a = 0.998. The central density of the torus p c = 10 11 cm" 3 . The 
ratio of the gas pressure to the total pressure J3 = 1.235 x 10" 5 . 



Differentiating the polytropic equation of state yields 
d a P = KTp T - l (d a p). 

The density structure of the torus is then given by 



(p 2 ' r p 

o a p = -a J — — + - — - 

H \ kT r - i 



(54) 



(55) 



Introducing a new variable £ (akin to logarithm of temperature), 
where £ = ln(T - 1 + Dep r-1 ), equation d55t simplifies to 



(56) 



The stationary and axisymmetry conditions imply that there are 
only two non-trivial components (r and 6) in the equation. By 
evaluating the line integral from r = r K , p = p c at the torus 



centre to the required (r, 6) location, the density field p(r, 6) is 
determined. 

If the pressure in the torus is dominated by the radiative pres- 
sure, we may set (1 - J3) ~ 1. Then we have 



P = fie 
kT = he 



45 



7i 2 {fim H pY 
45 



1/3 



n 4/3 



1/3 



,V3 



(57) 
(58) 



With p(r, 6) determined, the pressure and temperature are readily 
calculated. 

Figure 2 shows the density, pressure and temperature struc- 
tures of a model torus with n = 0.21, r k = 12r g ,/3 = 1.235xl0" 5 , 
and p c = 10 11 cm -3 , where r g = GM/c 2 is the gravitational ra- 
dius of a black hole with a mass M, and G is the gravitational 
constant. The spin parameter of the black hole a = 0.998 There 
are several noticeable features in the torus. The torus is pressure 
supported. Its boundary is located where the density p and hence 
the pressure P vanish. The temperature T also vanishes at the 
torus boundary surface for the equation of state that we adopted. 
As shown, the temperature T ~ 10 7 K in most of the torus in- 
terior, but it drops rapidly within a very short length of ~ 1 R g 
and reaches K at the torus boundary. Compared with the rota- 
tionally supported tori (Figure 1), the rotationally supported tori 
are larger in vertical extent. The vertical inflation of the torus 
is caused by the addition of radiative and thermal gas pressure 
forces. 

4. Radiative transfer calculations for accretion tori 

4.1. Emission from opaque rotationally supported tori 

For the opaque tori, the emission spectrum can be calculated 
from the emissivity distribution on the torus' boundary surface, 
with corrections for the relativistic shifts with respect to the dis- 
tant observer. We show in Figure 3 images of rotationally sup- 
ported tori viewed at inclination angles of 45° and 90° for vari- 
ous black hole spin parameters. The torus is left-right symmetric 
in shape if the black hole is not rotating. The black hole's ro- 
tation drags the surrounding space-time around the black hole, 
and so the torus around a Kerr black hole no longer has a left- 
right symmetrical shape. Optically thick tori would suffer self- 
eclipsing at high viewing inclination angles. Only the unob- 
scured regions on the torus surface would contribute to the emis- 
sion detectable by the distant observer. Each pixel in the torus 
images shown in Figure 3 has only one single value for the rel- 
ative frequency shift between the emission surface element and 
the observer. We coded this relative frequency shift with colours 
in the torus images. We can easily see that the regions with the 
highest redshifts and with the highest blueshifts are obscured by 
the front limb of the torus when the viewing inclination angle is 
close to 90°. 

In calculating of the intensity and the emission spectrum 
of the opaque tori, we employed a standard ray-tracing for- 
mulation similar to that used in the calculation of relativistic 
lines from geometr ically thin optically thick accretion disks (e.g. 
ICunninghamll975l) . since it adequately takes into account effects 
such as gravitational redshift, lensing, kinetic time dilation and 
Doppler boosting. Since only the propagation of radiation out- 
side the torus is relevant, we set the emission and absorption co- 
efficients to zero in the radiative transfer equation along the rays 
emerging from the torus' boundary surface. We computed the 
relativistic frequency shifts at the surface boundary, convolved 
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Fig. 3. False-colour frequency shift maps of the surface emis- 
sion from opaque tori around different black holes. The torus 
parameters are n = 0.232 and = 12 r g , where r g is the gravita- 
tional radius. The black-hole spin parameters are a = 0, 0.5 and 
0.998 (panels from top to botom). The viewing inclinations of 
the tori are 45° (left column) and 85° (right column). For a = 
the range of frequency shifts E/Eo for i = 45° and i = 85° 
are (0.874,1.445) and (0.756,1.560), respectively . Similarly, 
for a = 0.5 the frequency shift ranges are (0.870, 1.487) and 
(0.755, 1.591). Finally, for a=0.998 the corresponding frequency 
shift ranges are (0.864, 1.535) and (0.767, 1.616). 



this with a specified spatial profile for the source function of the 
emission and obtained the emission spectrum. 

Figure 4 shows the profiles of emission lines from an opaque 
torus and a geometrically thin, optically thick disk with the same 
inner and outer radius. The inner radius of the model accretion 
disk is much larger than the minimum radius allowed, which is 
about 1 r g for a maximally rotating black hole. A major differ- 
ence between a torus and a thin disk is that self-obscuration can 
occur for the torus at high viewing inclinations, while the first- 
order emission from the upper disk surface is always visible by a 
distant observer. Thus, for an opaque torus, self-eclipsing blocks 
the emission from the inner torus regions where relativistic ef- 
fects are most severe and the emission suffers the highest and 
lowest frequency shifts. For a thin disk, the inner disk regions 
with the highest and lowest frequency shifts contribute to the to- 
tal emission spectra. The line profiles of the torus and the disk 
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Fig. 4. Profiles of emission lines from an opaque rotationally 
supported accretion torus (top panel) and a geometrically thin 
optically thick accretion disk (bottom panel) viewed at different 
inclinations i. The torus dynamical parameters are n = 0.21 and 
rk = 12 r g . The inner boundary radius of the torus r in = 8.486 r g 
and the outer boundary radius r out = 20.246 r g . The disk has the 
same values for the inner boundary radius and the outer bound- 
ary radius as the torus. In both cases, the black-hole spin param- 
eter a = 0.998. The line emissivity is proportional to r~ 2 , where 
r is the radial distance from the central black hole. The line pro- 
files are normalised such that the flux F(Eq) = 1 at the view 
inclination angle i = 60°, where Eq is the rest- frame line centre 
energy. 



in Figure 4 show little difference at low viewing inclinations. 
This is easily understood, because the torus and disk not only 
look similar when they are viewed pole-on, they show similar 
radial dependences in the surface emissivity distribution and the 
entire upper emission surface obscured. For viewing inclinations 
i ~ 80° or higher, the line from the torus is narrower than the line 
from the disk, and in particular the edge of its blue wing is at a 
lower energy than that of the line from the disk, because the in- 
ner torus region where the bluest emission originates is eclipsed 
by the front limb. 

AAA. Emission from translucent pressure-supported tori 

Optically thin and translucent tori do not have a sharply defined 
emission surface. All parts within the tori contribute to the emis- 
sion, with the contribution weighted according to matter concen- 
tration (density) and to the local values of the thermodynamic 
variables (e.g. temperature) relevant to the radiation processes 
considered. The low optical depth across the torus permits the 
transmission of the emission from high-order lensed images, 
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therefore the emission spectrum of the torus is the sum of the 
spectra from all the orders of lensed images weighted by the op- 
tical depth of the ray. Figure 5 shows the intensity images of 
an optically thin radiative pressure dominated torus for different 
viewing inclination angles, and Figure 6 shows the intensity im- 
ages of an optically thin radiative pressure dominated torus with 
a lower internal radiation pressure. 

For the/? = 1 .235 x 10" 5 optically thin torus model, the value 
of J3 was chosen so that the inner and outer radii of the torus 
very closely matched the inner and outer radii of the opaque tori 
and thin disks discussed in the previous section, allowing a com- 
parison of the resulting images and emission line profiles. The 
second model with J3 = 5 X 10" 5 , i.e. a smaller internal radi- 
ation pressure, illustrates the dependence of the torus size on 
J3. The dynamical parameters of the two tori are n = 0.21 and 
rk = 12 r g , and the spin parameter of the central black hole is 
a - 0.998. The emissivity takes the form j oc p, and the light-of- 
sight optical depth across the tori r <*c 1. As seen in the figures, 
changes in the pressure ratio parameter J3 alter the aspect ratio of 
the torus, which determines to which degree self-eclipsing would 
occur for a given viewing inclination angle. In spite of this, the 
general emission properties as shown in the images are qualita- 
tively similar, because of the low optical depth in the torus. In 
both tori, the intensity of the emission is strongest at the interior 
of the torus, where the density is high. The rotation of the torus 
causes frequency shifts in the emission and Doppler boosting of 
the emission's intensity. These relativistic effects are more obvi- 
ous for high viewing inclination angles. The emission from the 
approaching limb of the torus is amplified and appears to be sig- 
nificantly brighter than the emission from the receding limb of 
the torus. 

Figure 7 shows the profile of an emission line from a ra- 
diative pressure dominated optically thin torus viewed at differ- 
ent inclination angles (top panel). The lines are broad and have 
an asymmetric profile, characteristic of line emission from rel- 
ativistic accretion disks viewed at moderate inclinati on ang les, 
i * 45° - 70° (see e.g. lCunninghanil975l:lFabian et all2000D . In 
contrast to the lines from relativistic accretion disks, the profiles 
of lines from optically thin tori do not change significantly when 
the viewing inclination angle changes from 45° to 85° (cf. the 
line profiles in Figure 4). The asymmetric broad profiles of the 
lines from optically thick relativistic disks are due to the com- 
bination of a number of effects: Doppler boosting, Doppler ve- 
locity shift, gravitational time dilation and gravitational lensing. 
Emission from an optically thin torus does not depend on the 
projected area of an emission surface element, as is the case of 
the optically thick accretion disks or optically thick accretion 
tori. While the emission from an optically thin torus is modified 
by Doppler boosting, Doppler velocity shifting and gravitational 
time dilation, it is less affected by gravitational lensing and the 
area projection effect. The transparency of the torus to the emis- 
sion could complicate the process of using the emission lines to 
diagnose the dynamical properties of the torus. First of all, the 
emission of the high-order lensed images, which are unobscured, 
could severely contaminate the emission from the direct image. 
Secondly, two emission lines can cross-contaminate each other 
in the spectrum, and because of the large relativistic shifts one 
line could cause absorption of another line that originates from 
a different region within the torus. For an optically thin torus, 
multiple rays corresponding to different energy shifts can hit the 
same pixel in the image plane. This is forbidden for ray-tracing 
in the case of optically thick disks or tori, where each ray orig- 
inates from a single point on the emission boundary surface of 
the disk or the torus and each pixel in the disk or torus image 



corresponds to a uniquely defined relativistic energy shift. We 
show in Figure 7 (bottom panel) that two distinct lines can be 
blended easily and may be taken as a single emission line with a 
complex relativistic line profile. 




Fig. 5. Surface brightness images of optically thin radiative pres- 
sure dominated accretion tori viewed at inclination angles of 15°, 
30°, 45°, 60°, 75° and 89° (left to right, top to bottom). The torus 
parameters are n = 0.21, = I2r g and/3 = 1.235 x 10" 5 . The 
black-hole spin parameter a = 0.998. The brightness of each 
pixel represents the total intensity integrated over the entire spec- 
trum. The torus brightness is normalised such that the brightness 
of the brightest pixel in each image is the same. 



4.1 .2. Emission from quasi-opaque pressure-supported tori 

The structured torus has density and temperature stratifications 
(see Figure 2) that give rise to opacity variation across the torus 
and opacity gradients along the line-of- sight for different radia- 
tive processes. When the density is sufficiently high, the torus 
becomes opaque to radiation. Quasi-opaque structured tori have 
more complex emission properties than their optically thin coun- 
terparts and the idealised optically thick rotationally supported 
tori in the previous sections. 

We investigated the emission properties using our covariant 
radiative transfer formulation, conducting a full radiative trans- 
fer calculation for a model- structured torus. We considered two 
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Fig. 6. Surface brightness images of optically thin radiation pres- 
sure dominated tori viewed at inclination angles of 30°, 45°, 75° 
and 85° (left to right, top to bottom). The torus parameters are 
n = 0.21, r k = I2r g and J3 = 5 x 10 -5 . The black-hole spin 
parameter a = 0.998. The brightness of each pixel represents 
the total intensity integrated over the entire spectrum. The torus 
brightness is normalised such that the brightness of the brightest 
pixel in each image is the same. 



opacity sources, whose specific emissivities in the rest frame are 
given by 



-2.5 



(59) 
(60) 



kT is the relativistic temperature and % 

-U™-3 u ,-1 



X 



where 

10" 46 erg s _1 cm" j Hz _i . C is a normalisation, dependent on 
P, chosen such that both 7*0,1 and 7*0,2 are equal at Eq = 0.1 
keV, at the torus centre. For p = 1.235 x 10 -5 , C = 2.162 x 
10" 45 erg s^cm^Hz -1 , whereas for£ = 5 x 10" 5 , C = 2.681 x 
10" 45 erg s _1 cm" 3 Hz _1 . The electron number density is defined 
as n e - plirniH- The corresponding specific absorption coeffi- 
cients are 



<*o,i(£o) = B i 



^0,2(^0) = B 2 



cm" 



cr T fi(E ) cm 1 



cm 



-3 



crrfiiEo) cm 



-1 



(61) 



(62) 



where crj is the Thomson cross-section, and f\(Eo) and /2CE0) 
are functions of photon energy. The torus is optically thin to the 
first process, but is partially opaque to the second process. We 
therefore set B\ = 0. Without loss of generality we consider an 
energy independent absorption, which implies fziEo) = 1. The 
normalisation constant B2 is chosen such that ao^r out -1-5 
across the torus, where r out is the outer boundary radius of the 
torus. Note that the first process has a similar density and temper- 
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Fig. 7. (Top panel) Profiles of an emission line from an optically 
thin pressure supported structured accretion torus viewed at dif- 
ferent inclinations i. The dynamical parameters of the torus are 
n = 0.21 and = 12 r g . The black-hole spin parameter a - 
0.998. The ratio of the gas thermal pressure to the total pressure 
P - 1.235 x 10" 5 . The radius of the inner boundary r in = 8.528r^ 
and the radius of the outer boundary r out = 20.246 r g . The line 
emissivity is proportional to the density p. These line profiles are 
normalised such that the line flux F(Eq) = 1 when the torus is 
viewed at i = 60°. (Bottom panel) Profiles of composite profiles 
from two emission lines. The torus and black hole parameters 
are the same as those for the lines in the top panel. The line en- 
ergies are such that one line has an energy 10% higher than the 
other line and the emissivity of the line with the higher line cen- 
tre energy is 14% of that of the line with the lower line centre 
energy (c f. analogous to the relative properties of the Fe Ka and 
Kp lines dHolzer et al.ll 19971) ). 



ature dependence to thermal free-free emission. The second pro- 
cess mimics a free-electron scattering-like process that converts 
photons with different energies indiscriminately into a power- 
law energy distribution. 

Figures 8 and 9 show intensity images of opaque radiative 
pressure dominated tori using the aforementioned spectral pa- 
rameters. These tori differ from the optically thin case in that the 
emissivity and opacity have changed, causing significant differ- 
ences in the intensity images. Towards the outer-edges of the tori 
there is now obvious limb darkening (dark red). There is now 
an additional temperature dependence on the emission that is 
caused by the temperature stratification of the torus, nearly tan- 
gential rays sample only the cooler surface layers of the torus. 
However, rays almost perpendicular to the torus surface have a 
much higher probability of traversing the hotter layers deeper 
beneath the torus surface. Consequently, surfaces viewed face- 



10 



Younsi, Wu and Fuerst: General relativistic radiative transfer 




Fig. 8. Surface brightness images of opaque radiative pressure 
dominated accretion tori viewed at inclination angles of 15°, 30°, 
45°, 60°, 75° and 85° (left to right, top to bottom). The torus pa- 
rameters are n = 0.21, r k = I2r g and J3 = 1.235 x 10" 5 . The 
black-hole spin parameters = 0.998. The emissivity is provided 
by emission from two spectral lines, and the two opacity sources 
given in equations (59) and (60). Absorption is provided through 
the Thomson cross-section in equation (62). The brightness of 
each pixel represents the total intensity integrated over the en- 
tire spectrum. The torus brightness is normalised such that the 
brightness of the brightest pixel in each image is the same. 



on by an observer will appear to be much brighter than those 
viewed at higher inclination angles. For these tori we have as- 
sumed r » 1 . 

4.2. Some remarks 

4.2.1 . Optically thin vs optically thick emission 

For the optically thin tori in figures 5 and 6, the low optical 
depth means that these tori are almost transparent, allowing the 
entire torus volume to contribute to the emission. The emission 
from the approaching side of the torus is much stronger due 
to Doppler boosting, which causes an increase in the pro- 
jected velocity of the gas along the line of sight towards the 
observer. The receding limb's projected velocity is decreased 
and consequently dimmed. At higher viewing inclination angles 
geodesies travel through more emissive material than at lower 




Fig. 9. Surface brightness images of quasi-opaque accretion tori 
around extreme Kerr (a=0.998) black holes, viewed from left to 
right, top to bottom, at observer inclination angles of 30°, 45°, 
75° and 85°. Same parameters as figure 8, but now J3 = 5 x 10" 5 . 
Line emission and continuum emission are included. These tori 
also self-occult and higher-order emission is greatly suppressed 
at high inclination angles. The images are normalised such that 
the brightest pixel in each image is of the same intensity. 



inclinations. Since the fast-moving gas in the inner region of the 
torus is visible, blueshifting and beaming are very strong and 
emission from the approaching side of the torus dominates these 
images. The overall brightness of the optically thin tori at high 
inclination angles is much brighter than at lower inclinations, 
because at high inclinations geodesies may sample significantly 
more emissive material in the limbs. 

For the rotationally supported opaque tori, emission comes 
from the surface of the torus only. Since the opaque torus 
models are optically thick, the front limb of the torus obscures 
the fast-moving inner regions of these tori at higher inclination 
angles. This is reflected in their emission spectra in figure 4, 
where for i > 45° the observed flux decreases rapidly. For the 
optically thin torus model, emission is observed from the entire 
volume of the torus, not just the surface. Since the fast-moving 
gas from the inner-region contributes to the observed emission 
line profile (Figure 7), this gives rise to the more pronounced 
red and blue wings, which are observed even at high inclination 
angles. Additionally, for optically thin tori the observed line 
profiles appear to be roughly monotonically increasing with 
observer inclination angle. This is in contrast to the line profiles 
from optically thick accretion disks and tori in Figure 4, which 
decrease in amplitu de and appear wedge-shaped for inclinatio n 
angles beyond 45° fl Kojimal (Il99lh . iBeckwith & Donel (l2004bh . 
iFuerst & Will (I2004 ). As was found in iFuerstl (l2006h . altering 
the black hole spin does not affect the observed torus images 
or emission line profiles significantly. The inner edges of the 
tori presented in this paper do not extend as close to the event 
horizon as thin disks. Consequently, the red wing of their line 
profiles is not as extended. The shape of these tori, as well as 
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the location of their inner edges, depends much more on the 
distribution of pressure forces within the tori themselves. If 
these models are realistic, it would prove difficult to derive the 
central black hole spin from the observed spectral information. 

The emission from quasi-opaque tori is restricted to a thin 
surface skin layer that is similar to cases of the opaque accretion 
tori and disks. However, the quasi-opaque torus is distinguish- 
able from the fully opaque accretion tori and accretion disks by 
the effects arising from gradients in the optical depth. This effect 
is particularly noticeable near the edges of the quasi-opaque 
tori, where significant limb darkening occurs. These limb effects 
are a manifestation of the temperature and density stratification 
within the tori. When looking at the edge of a torus, for a given 
optical depth, cooler, less dense material is observed than when 
looking towards the torus centre. Just above the surface of the 
quasi-opaque torus is an optically thin photosphere. As can be 
seen in figures 8 and 9, photons traversing only the edges of the 
limb of the torus are weakly attenuated by the cool material and 
so appear significantly darker than the rest of the torus surface, 
i.e. only line photons are observed. Hotter neighbouring regions 
appear to be much brighter because continuum photons are 
attenuated by the local emission and absorption properties of 
the emitting medium. 

We illustrate the limb effect with a simple phenomenolog- 
ical model of an emitting and absorbing medium, in local 
thermodynamic equilibrium (LTE), at a uniform temperature T, 
with a background source of intensity 1(0). In the absence of 
scattering, taking To = at the outer edge of the limb, equation 
(18) implies 



to energy redshifts, whereas photons propagating in the opposite 
direction, deeper into the gravitational well, are blueshifted. 



J(r v ) = J(0)e-^+^(r)(l-e- T 0, 



(63) 



where !B(T) is the Planck function. For low optical depths, 
r y — > and so J(r v ) — > 1(0), i.e. the observed intensity of the 
limb tends towards the background intensity, is independent of 
temperature, and emission from regions deeper within the torus 
becomes negligible. In our calculations we assumed 1(0) <?c 1 
and so the limb darkening effect is obvious, and in strong con- 
trast to the rest of the torus image. 



4.2.2. Gravitationally induced line resonance in 3D flows 

Consider a two-level atomic system. Absorption of a photon 
causes the excitation of an electron from a lower energy state 
to a higher energy state. Conversely, a photon is emitted when 
an electron is de-excited from a higher energy state to a lower 
energy state. The energy of the absorbed and the emitted photon 
is the same as the energy difference between the two states. Such 
a system can be considered an oscillator. Under general condi- 
tions, two oscillators with different intrinsic frequencies at dif- 
ferent locations would not be resonantly coupled, and two lines 
from different atomic transitions would not exhibit resonant be- 
haviour without coupling to either other lines or an optical pump. 

However, in the vicinity of black holes, where relativis- 
tic effects are severe and gravity is extreme, lines from differ- 
ent atomic transitions can couple and exhibit re sonance phe- 
nome n a. As shown in num e rous publicati o ns (e.g.lCunninghaml 
19751: lFabian etal.1 1 19891 IStellal Il990h iFanton et al.1 119971: 
Reynolds et al.l!999b Fabia n et al.l2000l:lFuerst & Wul2004 . the 



gravitational frequency shifts of lines from relativistic accretion 
disks around black holes can be severe. Photons propagating up- 
wards, out of the gravitational well of the black hole, are subject 



■ E 2 



T2 



Ti 



■E 2 ' 




Black Hole 



r 2 



ri 



Fig. 10. Schematic illustration of the resonance between two 
lines with E\ and E2, emitted at radial distances r\ and from 
a black hole, where r^> r\. In the rest frames E\ > E2. 



Consider a radiative transition process occurring in a 
medium at a radial distance r\ from the black hole, which con- 
sequently emits a photon of energy E\ . The photon propagates 
outward, and its energy is redshifted. When the photon reaches 
a distance ri, its energy becomes E2, lower than E\, the energy 
at its point of emission. Suppose that the photon encounters an 
electron and is absorbed. The electron is excited to a higher en- 
ergy state. The electron is then de-excited and returns to its orig- 
inal state, emitting a photon of energy E2. This photon, how- 
ever, propagates inward. When it reaches r\ , its energy has been 
gravitationally blueshifted and is now E\. Again suppose that 
this photon encounters an electron, is absorbed, and causes the 
electron to be excited to a higher energy level. This electron is 
subsequently de-excited, and again emits a photon of energy E\ . 
The photon propagates outward once more, is absorbed by an 
electron, and causes another excitation. The later de-excitation 
of this electron leads to the emission of another photon of en- 
ergy £2, and the photon propagates inward again... This process 
would persist, forming a resonant feedback cycle between the 
transitions of the two lines with different rest- frame energies (see 
Figure 10). 

Such a line resonance occurs easily in 3D flows but not in 
geometrically thin accretion disks or 2D flows. This can be un- 
derstood as follows. In the Schwarzschild space time, for a pho- 
ton with an energy E\ located at r\ , one can also find a closed 
surface corresponding to an energy redshift of AE, such that the 
photon energy becomes E2 - E\- AE, located at r2 > n . The 
converse also holds, and a photon of energy E2 located at r2 al- 
ways finds a closed surface corresponding to a blueshift of AE, 
such that the photon energy increases to E\ - E2 + AE. A pho- 
ton emitted from the vicinity of a black hole must pass through 
the closed surface with a specific energy shift before reaching 
infinity, and in 3D radial flows the entire surface is in principle 
embedded completely in the flow. A 2D flow cannot completely 
contain such a surface and hence the photon can escape to infin- 
ity, passing through the surface with a specific energy shift at a 
location not inside the flow. The situation is similar in the Kerr 
space-time. A detailed quantitative analysis of such a line res- 
onance, and relativistic radiative transfer calculations including 
these masing effects warrants further investigation. We leave this 
to a future article. 
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5. Conclusions 

We have derived the radiative transfer equation from first prin- 
ciples, conserving both particle number and phase- space den- 
sity. The equation is thus manifestly covariant, and it is applica- 
ble in arbitrary 3D geometrical settings and for any pre-defined 
spacetime metric. The more general form of the equation ex- 
plicitly considers a covariant particle flux with mass. We have 
found that ind addition to modifying the geodesic trajectories 
that connect the observer and the emitting regions, the presence 
of particle mass in the radiative transfer equation introduces a 
mass-dependent aberration, which reduces the intensity gradient 
along the ray. In the zero-mass limit this general radiative trans- 
fer equation recovers its original form for the massless particles. 

We carried out demonstrative numerical general relativistic 
radiative transfer calculations with a ray-tracing algorithm con- 
structed from the formulation. Different 3D accretion tori around 
rotating black holes with different geometrical aspects, physical 
structures, emission properties, and optical depth variations were 
considered. We demonstrated that radiative transfer calculations 
based on the formulation are able to deal with the complexity 
in the various combinations and convolution of relativistic, ge- 
ometrical, physical and optical effects. Our calculations clearly 
showed the significant role that structures and optical depth, and 
their gradients, together with geometrical and relativistic factors, 
play in shaping the emission properties of these 3D relativistic 
flows in the vicinity of rotating black holes. The calculations also 
showed the presence of limb effects in 3D objects with finite op- 
tical depths. 

We note that gravitationally induced line resonance can oc- 
cur in 3D accretion onto a compact object. This phenomenon is 
not present in 2D planar objects, such as geometrically thin ac- 
cretion disks, where the radiation can escape from the disk sur- 
face to free space without additional absorption or re-emission. 
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